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We analyze the consequences of Lorentz violation (LV) to three-generation neutrino 
oscillation in the massless neutrino sector. We present a general formalism of three-family 
neutrino oscillation with neutrino flavor states being mixing states of energy eigenstates. 
It is also found that the mixing parts could strongly depend on neutrino energy by special 
choices of Lorentz violation parameters. By confronting with the existing experimental 
data on neutrino oscillation, the upper bounds on LV parameters are derived. Because 
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energy dependence may test and constrain the Lorentz violation scenario for neutrino 
oscillation. 
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1. Introduction 

Lorentz invariance has been a fundamental principle in physics since Einstein estab- 
lished special relativity in 1905. In the last century, numerous experiments provided 
precise verifications of Lorentz invariance, and most current data are consistent with 
this symmetry. Over the last decade there is a growing interest in studying Lorentz 
invariance violation in the physics society. Though no theoretical model predicts 
Lorentz violation conclusively, there have been some theoretical suggestions that 
Lorentz invariance may not be an exact symmetry at all energies. In particular, 
many works to describe the force of gravity within the context of a qu antum theory 
imply the breaking of Lorentz symmetry (including string theory ^1^, warped brane 
worlds'^, and loop quantum gravitj^)^. Other high energy models of space structure 
could also contain Lorentz violation^. Even if Lorentz symmetry is broken at high 
energy, there can still be an attractive infrared fixed point . So we can still get an 
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approximate Lorentz invariant world at low energy range. Besides above theoret- 
ical motivations, many ot her models, such as emergent gaugebosonJ^E!^ varying 
moduli^, ghost c ondens atj^, space-time varying couplin j^^ l ^^ [ or varying speed 
of light cosmologjff^S^l also incorporate Lorentz violation. Whether Lorentz sym- 
metry is perfectly unbroken under all conditions is still an important theoretical 
question. 

Lorentz invariance violation was proposed as a solution to two important ex- 
perimental problems, i.e., the observation of TeV photons and of cosmic ray events 
above the GZK cutoffl^ ^'^ , for a review see Ref. HDl In the first case, it 
is hard to observed ultrahigh energy photons coming from ulterior galaxies since 
cosmic gamma rays with energy above 10 TeV should interact with cosmic infrared 
background photons and convert into electron-position pairs. However, the High 
Energy Gamma Ra y A stronomy (HEGRA) satellite detected 24 TeV gamma rays 
from Markarian 5oJ^. In the second case, ultra high energy cosmic rays interact 
with cosmic microwave background photons and produce pions. The cosmic rays 
lo se ener gy through this process until they pass below the GZK cutoff energy 10^^ 
e\ 122 | 23 | However, there has been report that the cosmic ray spectrum extends be- 
yond this energy24|_ Lorentz violation could give a preparative solution to the two 
problems because the threshold energy at which the cutoff occurs could be altered 
by modifying special relativity, though it is n ot the only way t o ex plain these prob- 
lems. We also awarJ^ that recent HiReJ^ and Pierre AugeJ^ measurements of 
the ultrahigh energy cosmic rays show a sharp suppression around the energy of 
10^^ eV, which is consistent with the expected cutoff energy. If the phenomenon of 
GZK cutoff is confirmed by further experiments, Lorentz invariance will still be a 
good symmetry in the area of ultra high energy cosmic ray. Therefore one would 
need to look at other possibilities for the breakdown of Lorentz invariance. 

Neutrinos provide another interesting laboratory for studying the possibility of 
Lorentz violation. If Lorentz invariance is violated by quantum gravity, the natu- 
ral scale one would expect a strong Lorentz violation is the Planck energy of 10^^ 
Qevfm, But the attainable energy of accelerator is of TeV scale. So it is difficult 
to directly detect Planck scale Lorentz violation in the laboratory. However, there 
might be a small amount of Lorentz violation at lower energies if Lorentz sym- 
metry is violated at Planck scale. Neutrinos offer a promising possibility to study 
Lorent z vi olation that may exist at the low-energy as the remnants of Planck-scale 
Physic*^. Neutrino oscillation is an important problem in neutrino physics. To 
explain this kind of phenomena, the conventional scenario is to assume that neutri- 
nos have masses. In this assumption, there is a spectrum of three or more ne utrino 
mass eigenstates and the flavor state is the mixing state of mass eigenstateJ^^^^. 
Coleman and Glashow pointed out that neutrino oscillation can take place even 
for massless neutrinos if Lorentz invariance is violated in the neutrino sectoJ^. In 
Ref. 16, they assumed that the maximum attainable speed of a particle depends 
on its identity and that the flavor states are the mixing states of speed eigenstates, 
and analyzed neutrino oscillation quantitatively. In Ref. |33| a two-generation model 
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to study neutrino oscillations is established. This model involves a mass term and 
a single nonzero coefficient for Lorentz violation. Kostelecky and Mewes presented 
a general formalism for violations of Lorentz and CPT symmetry in three-family 
massive neutrino oscillatiori24l. in their calculation the mixing part has the same 
form comparing with that of the conventional massive neutrino scenario because of 
their assumption that neutrino are massive. They also built a special model called 
bicycle model in massless neutrino sectoJ^, in which the Lorentz violation param- 
eters are direction-dependent. Ref. |36] built a three-neutrino massive model with 
Lorentz-violating terms. All classes of neutrino data are descried in this model, in- 
cluding LSND oscillation data. Grossman et al. studied the interactions between the 
neutrino and the Goldstone boson of spontaneous Lorentz violation, and proposed a 
novel dynamic effect of Goldstone-Cerekov radiation, where neutrinos mo vin g with 
respect to a preferred rest frame can spontaneously emit Goldstone bosonJ^Il. Arias 
et al. analyzed the consequences of Lorentz violation in the massless ne utri no sector 
by deforming the canonical anti-communication relations for the fields Morgan 
et al. analyzed atmospheric neutrino oscillations at high energy by modified disper- 
sion rela tion s and placed bounds on the magnitude of this type of Lorentz invariance 
violatioi|22l. 

In this paper, we study Lorentz violation contribution to neutrino oscillation. In 
our calculation we assume that neutrinos are massless and that the neutrino flavor 
states are mixing states of energy eigenstates. We calculate neutrino oscillation 
probabilities by the effective theory for Lorentz violation, which is usually called the 
standard model extension (SME)^*^ . In our work, the mixing angles for neutrinos 
are functions of Lorentz violation parameters. 

The structure of this paper is as follows. In Sect. [51 we figure out the neutrino 
oscillation probabilities by the effective theory for Lorentz violation. In Sect. [31 we 
introduce some specific models and give the numerical values for LV parameters by 
comparing our theoretical oscillation probabilities with experimental results. Then 
we analyze the dependence of the new oscillation equation on neutrino energy and 
neutrino propagation length. Remarks and conclusions are given in Sect.[4l 

2. Theoretical framework for Lorentz violation and neutrino 
oscillation probabilities 

In the conventional massive neutrinos scenario, neutrino mass eigenstates are com- 
ponents of neutrino flavor states. Neutrinos change from one flavor to another dur- 
ing the propagation, because the phase of time for each mass eigenstate is different. 
However, Lorentz invariance violation may be another origin for neutrino oscilla- 
tion. In this work, we analyze the consequences of Lorentz violation in the massless 
neutrino sector. Different from the conventional massive neutrino model, neutrino 
flavor states are mixing states of eigenenergy in our calculation. The bounds on 
LV parameters will be given by comparing our calculation with the experimen- 
tal results. In this part we will flgure out the neutrino oscillation probability. The 
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general and detailed calculation method for neutrino oscillation was proposed in 
Ref. mi However, for the integrity of the article, we display the entire calculation 
for neutrino oscillation probabilities. ^^^^ 
In the framework of standard model extension, we consider the LagrangianlSlIl2l 
for neutrino sector given by 

c = \f.i^dlyj,^ + \ic^:^y.i''y^vu ~ ^sT^A^^s, (1) 

where /i, v are spinor indices and A, B are flavor indices, c^^ and a^^ are LV 
parameters. In Eq. ([1]), the first term is consistent with the minimal standard model; 
the second and third terms describe the contribution from Lorentz violation, which 
denote CPT even term and CPT odd term respectively. For the convenience of 
calculation, we transform the lagrangian as 

where C equals to L in quantum field theory, because divergence terms have no 
contribution to the action. The Euler-Lagrange equation of motion for neutrinos 
can be written as 

z7°5o^^, + ifd^v, + icZl'^d-v, - a^:^7^^« - 0. (3) 

A new motion equation can be obtained by multiplying Eq. ([3]) with the matrix 7*^. 
Comparing the new equation with the conventional equations of motion {iS^sdo — 
T-Lj^b)i^a = 0, we obtain the Hamiltonian for neutrinos 

n = ~i-i°fd, - ic^^7°7''a'^ + a'j^7°7^. (4) 

Note the unconventional time-derivative term in Eq. ^ has been included in Eq. 
The above discussion is within the context of quantum field theory. Now we trans- 
form the description into quantum mechanics and treat the differential operator as 
momentum operator to study neutrinos with a fixed momentum. In this paper, we 
focus on the upper bounds of Lorentz violation contribution to neutrino oscillation 
and so far right-handed neutrinos or left-handed antineutrinos have not been de- 
tected experimentally. So we assume that neutrinos are massless. Our calculation 
is based on three-generation model. With the basic vector (it^(p), p))-^, the 
Hamiltonian matrix for neutrinos can be given by 





'■^^ \ P\ --"IPI 

Substituting Eq. ([5]) into Eq. ([3|) the general dynamical equation for neutrinos 
could be obtained. Then we can study the dynamical character of neutrinos by 
calculating the dynamical equation. However, we are more interested in neutrino 
oscillation probability in this work. So we just need to figure out the eigenenergy 
for neutrinos by diagonalizing the Hamiltonian. In Eq. the up-diagonal term 
is the Hamiltonian for left-handed neutrinos and the down-diagonal term is the 
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Hamiltonian for right-handed antineutrinos. We just study the left-handed neutrino 
part because the two terms have the similar form. We simplify the Hamiltonian -ff^^ 
as 

= (6) 

where B — e, /i, r. To figure out neutrino eigenenergy we diagonalize the Hamil- 
tonian matrix ([6]) by a 3 x 3 matrix U . Note that [/ is a unitary matrix. The 
eigenenergy matrix is 

E,,=u]^h,,U,,, (7) 

where is a diagonalized matrix. The eigenenergy can be labeled as , where / = 
1, 2, 3. Neutrino energy eigenstates are the linear combination of flavor eigenstates 
because matrix U is unitary: 

where / and A represent different eigenstates and flavor eigenstates respectively. 
Because the matrix U is unitary, Eq. ([5]) can be transformed as 

\^At)) ^Y.^u')]^e-^^^\u'),,\v,). (9) 

IB 

Neutrinos propagate at the speed of light approximately because we have assumed 
that neutrinos are massless and that LV parameters is small. When the propagation 
length is L, the flavor state can be written as 

\vAL)) = Y.^u')]^e-^^'\u'),,\v,). (10) 

IB 

^From Eq. (fTO|l and the unitary of U , the oscillation probability can be expressed 

as 

Pi^A -^Vb) = ^AB 

-4^5R[([/^):j{/\j{/V^^(C/^,,]sin2(^L) 

i>j 

+ 2^3[([/^);j{/"^),Jt/^)*Jt/^),,]sin2(Ai?,,L), (11) 
i>j 

where 3? and denote the real and imaginary parts respectively. 
3. Specific models and bounds on LV parameters 

In Sect. [51 we obtained the analytic equations for neutrino oscillation. In this section, 
we discuss three specific models. Comparing with the experimental results, we will 
figure out the upper bounds on LV parameters in different models. The Lorentz non- 
invariant direction-depe ndent oscillations for massless neutrinos are not supported 
by recent research wor]J22I311_ go in this part we try to build a massless neutrino 
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model without direction-dependent oscillations. All the three special models in our 
work are direction- independent. 

In Eq. ([6]), is a 3 x 3 matrix in neutrino generation space. To reduce the LV 
parameters, we assume that the neutrino Hamiltonian can be simplified as 

/E e 0\ 

h,,^ i e E + 7J C] , (12) 

\o C eJ 

where i? = |^|, e = ^'^.p^+a^ = c'^^E + a^^^, C = c'l^^E + a"^,, and 77 = c^^^E + a';^^. 
Lorentz violation parameters are defined in a sun-centered inertial frame. There are 
6 non-zero LV parameters in the above specific model, and the other LV parameters 
are zero. 

By substituting Eq. (I12p into Eq. ([7]), the eigenenergy matrix E and mixing 
matrix U could be figured out as 



/E 







E,,^ E+- 







""^h- — 2 — -I 



(13) 



and 



where 



c 



u = 



e 
N 



ry- v/4(£2 -I- (2) + 772 ( 



2VN 



N 



2^/M 



M 



(14) 



2 , .2N , , f?v/4(e2 + C2)+77' 



M = 2(^2 + ^2) + ^ + 



2 , /■2^ , 77^4(£2+(2)_^^2 



N = 2(e^ + C') + 



(15) 



With the eigenenergy matrix E and the unitary matrix U, the oscillation probabil- 
ities can be expressed as 
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P{Pe ^ i^e) = 1 



■ sm [(- - )L\ 



sin [(- + ^ )L] 



4£4 . ^ ^4(£^ + C^)+??^ 



, 4£^ . 2,(V^^^TW+^,., 
^(-e - = ^^T^^^^sm [( )L], 



P(;.^ ^ ^.^) = 1 _ -^—^^—^sin [( ^ )L], 



T>f ^ 4C^ . 2r, V4£2 + 4C2+r;2 



p. . 4£V . n ^4^2+4^2+^ 

P(l/e Z^r) = , sm [(- - -^^ )L] 

4£V .sin^[(^ + ^i^!±i^)L] 



(£2 + C^)[v/2(£2 + C^) + f + 



(£2 + C2)(4£2+4C2+772) 



4£2C2 . ,^^ V4£2+4C2+y?2 

sm [( )£], 



4£2C2 . 2 77 v/4£2 + 4C2 + 772 

■sm [(7 z )L] 



(£2 + C^)[V2(£2+C2) + f - "^^'^y^'^^'^ 



4£2C2 . 2r.^ , V4£2 + 4C2+r72 

■ sm [(7 + )L] 



(£2 + e)[^2is^ + e) + 'i + '^"^'Y^'''^ 



s^^li^^W^n- (16) 



(£2 + C2)(4£2+4C2+772) 2 

In Eq. (|16p . oscillation probabilities satisfy the probability unity equation 
"YligPi^A ^ ^b) — 1; which is guaranteed by the unitarity of U. From above 
calculation, we can see that not only the eigenenergy splitting A E^j but also the 
mixing parts are functions of LV parameters. To clearly understand the novel effect 
of Eq. Uni), we can analyze P{v^ Vt) with a special choice of LV parameters. If 77 
is the only nonzero CPT even LV parameter, then 7/ = -E, £ — aP^^^ C = oP^t- Now 
both the mixing part and the eigenenergy splitting A Eij are functions of neutrino 
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energy. In the conventional massive neutrino model, the corresponding probability 
can be expressed as P(t'^ — > i't) = sin^ 26* • sin^[1.27 Am^(L/£^)], while the mix- 
ing part sin^ 26 is independent of energy. So the amplitude of conventional massive 
neutrino oscillation is fixed, and the phenomenon of oscillation disappears when the 
neutrino energy is high enough because the oscillation P directs to L/E at high 
energy. But in Eq. ([16]), we can clearly see that the mixing part of Vr) is a 

function of neutrino energy. So oscillation amplitude varies with neutrino energy. By 
further calculation, as shown in l3.1l the oscillation probability will be suppressed at 
high energy, which is consistent with the conventional massive model, but the sup- 
pressed threshold energy is much higher comparing with the conventional massive 
model. 

There would be different novel behaviors of the neutrino oscillation in the pres- 
ence of Lorentz violation by choosing different LV parameters. In the following we 
will study three special cases with different kinds of choosing LV parameters. By 
comparing with the exiting experimental data, we can constrain LV parameters. 

3.1. Model 1 Non-diagonal terms are CPT-Odd parameters with 
CPT-Even parameters only in the diagonal term 

If we assume that Og^, aj^^, and c^^^^E are nonzero, then 

£ = 0°^, C-a°., ^ = (17) 

Substitute these LV parameters into Eq. pTjl . To calculate the numerical results of 
LV parameter, we compare the theoretical equations with experimental results. 

KamLAND detected P{ue —> Ve) — 61% with neutrino energy E ~ 4.3 MeV 
and neutrino propagation length L ~ 180 kmS^. MINOS observed muon neutrino 
disappearance with neutrino energy E ~ 4.9 GeV and propagation length L ~ 735 
km, P{Vfi v^) ~ 76% We also get the flavor change channel of muon neutrino 
to tau neutrino from K2K, P{v^ Vr) ^ 36% with E ~ 1.8 GeV, L ~ 250 kmE3. 

Substituting these experimental data into Eq. (|lip . we get three equations with 
three LV parameters. Because this is a nonlinear system of equations, the solution 
is not unique. But we can give the order of LV parameters in principle: a^^ ^ aP^^ ^ 
10-" eV and 0°°. ~ 10-2°. 

For further restriction of LV parameters, we use the experimental results for 
the oscillation channel of electron neutrino to muon neutrino. LSND has used muon 
sources from the decay 7r+ fJ-'^ to detect neutr ino oscillation in the subsequent 
decay of the muon through /i+ e+ + -I- i/^ This experiment finds the 
oscillation -> Ve (with 20 MeV < E,,^ < 58.2 MeV) with a probability of 0.26%. 
However, the MiniBoone experiment, which was expected to verify the results of 
LSND, has reported their first result which does not favor the simple explanation 
of the LSND results based on the two flavor neutrino oscillation '^^ . In this paper, 
we use the experimental results from K2K about the oscillation channel i/^ — )■ v^. 
with the mixing angle sin 26i^e < 0.13 Now we plot the oscillation probability 
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as a function of the neutrino energy i?, and as a function of the path length L, 
respectively, with the LV parameters a^^ = 1.09 x 10"-^^ eV, = 2.97 x 10~^^ 
eV, and c™ ^ 1.42 x lO^^o. 



? 0, 




> 0.8 - 



100 200 300 

E/IOMeV 



(a) 



100 

L/km 
(b) 



Fig. 1. Neutrino oscillation probabilities for — » i^e- (a) The left plot shows the oscillation 
probability as a function of neutrino energy with a fixed path length L = 180 km. (b) The right 
plot shows the oscillation probability as a function of path length with a fixed energy E = 10 
MeV. 




Fig. 2. Neutrino oscillation probabilities for — > u^. (a) The left plot shows the oscillation 
probability as a function of neutrino energy with a fixed path length L = 100 km. (b) The right 
plot shows the oscillation probability as a function of path length with a fixed energy E = I GeV. 



Fig. [T] shows the oscillation probability for i/e — > I'e as a function of the neutrino 
energy E with a fixed path length L = 180 km, and as a function of the propagation 
length L with neutrino energy E = 10 MeV, respectively. Comparing with the re- 
sults of massive neutrino scenario, we find that in Lorentz violation model electron 
neutrino does not have drastic oscillation at low energy. This is different from the 
conventional massive neutrino scenario. In conventional neutrino model, the phe- 
nomena of neutrino oscillation (j/g — > j/g) disappears at high energy {E > 100 MeV 
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Fig. 3. Neutrino oscillation probabilities for U/j, — > 1/^. (a) The left plot shows the oscillation 
probability as a function of neutrino energy with a fixed path length L = 100 km. (b) The right 
plot shows the oscillation probability as a function of path length with a fixed energy E = I GeV. 



with a fixed path length L = 180 km). However, in Fig. [T] we find the existence of 
oscinations at high energy. This distinction is resulted from different dependence on 
neutrino energy for oscillation probability. In addition, the mixing part is fixed in 
convention massive model. But in Fig. [5] we could find that the neutrino oscillation 
amplitude varies with neutrino energy. These novel effects can be tested at high 
energy sector. Because the dependence on path length has the same form in two 
different models (massive neutrino model and Lorentz violation model), the plots 
of path length dependence have no remarkable distinction with the conventional 
massive neutrino model. 



3.2. Model 2 non-diagonal terms are CPT-Even parameters with 
CPT-Odd parameters only in the diagonal term 

If we assume that aj^^, CgJ^, and cj^*^ are nonzero, then 

e = c™i?, C = c°°i?, v = al^. (18) 

Substitute these LV parameters into Eq. (fTTj) . Similarly to Model 1, we use 
KamLAND, K2K and MINOS experimental results to compute the parameters. The 
numerical values are also consistent with experimental results about the oscillation 
channel i/^ — > ly^ from K2K. Fig. [3] shows the oscillation probability of the channel 
j/^ — !■ t/^ as a function of the neutrino energy E and as a function of the propagation 
length L respectively with LV parameters a"^ = 1.58 x 10~^^ eV, cf^ ~ 2.93 x 10~^^, 
and cO° = 7.68 x lO^^^. 

The order of LV parameters c^^ is larger than that of Model 1, so the oscilla- 
tion periods of the changing probabilities against neutrino energy and path length 
are much smaller than those in Model 1. In addition, the oscillation amplitude is 
steady at high energy, which is different comparing with Model 1. When neutrino 
energy is high enough, parameters s and C are much lager than rj because e and ^ 



6, 2009 13:15 WSPC/INSTRUCTION FILE neutrino 



Lorentz violation in three-family neutrino oscillation 11 

linearly depend on neutrino energy. LV parameters c™ and cj^° in the numerator 
and denominator have the same power in the mixing part, so the values of mixing 
part are steady at high energy. 

3.3. Model 3 Both of CPT-Even and CPT-Odd parameters are 
included in the diagonal term 

If we assume that a^^, a°^, aj^^, and cj^^ are nonzero, then 

£ = C = «°r, '7 = < + c™i?. (19) 

We introduce a new equation C = a;e, then the constraint of the mixing part 
sin^ 20pe < 0.13 from K2K experimental results can be easily satisfied by adjusting 
the constant x. In our calculation we choose x — 2.7. Substitute all these new pa- 
rameters to Eq. (fTT|) . Comparing with KamLAND, K2K and MINOS experimental 
results, we can figure out the LV parameters. Figs. [4] and [5] show the oscillation prob- 
ability against neutrino energy E and propagation length L respectively with LV 
parameters 0°^ = 1.28 x 10~" eV, a" = 3.40 x 10"" cV, and c°" = 1.04 x lO^^o. 




0,5 I . 1 . 1 ^ 1 . 1 . 1 0.5 I ■ ' ■ ' ■ ' ■ 1 

100 200 300 400 500 50 100 150 200 

E/10MeV L/km 
(a) (b) 

Fig. 4. Neutrino oscillation probabilities for — » Ve- (a) The left plot shows the oscillation 
probability as a function of neutrino energy with a fixed path length L = 180 km. (b) The right 
plot shows the oscillation probability as a function of path length with a fixed energy E = 10 
MeV. 

Figs, m and [5] keep the main characteristics of Figs. [T] and [2l Neutrinos do not 
have drastic oscillation at low energy and the oscillation phenomena still exists at 
high energy, which is different from the massive neutrino scenario. Similar to Model 
1, the oscillation amplitude varies in different energy scale. These similarities can be 
explained by the fact that the oscillation probability dependence on LV parameters 
and neutrino energy have similar form in the two models . 

During the above calculation the boost effect has not been considered. In this 
work we assumed that the LV parameters are direction independent in the sun- 
centered frame. The experiments are done on the Earth which is moving around the 
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0.0 1 ' 1 ■ 1 ■ 1 ' 1 ■ 1 0.2 1 ' 1 ■ 1 ■ 1 ' 1 ■ 1 

2 4 6 8 10 20 40 60 80 100 

E/GeV L/km 
(a) (b) 

Fig. 5. Neutrino oscillation probabilities for j/^i — > i/^. (a) The left plot shows the oscillation 
probability as a function of neutrino energy with a fixed path length L = 100 km. (b) The right 
plot shows the oscillation probability as a function of path length with a fixed energy E = 1 GeV. 



sun and is also rotating. The experiments in laboratories therefore see LV parameters 
which have direction dependence because there is a boost. When the long baseline 
experiments based on the earth are discussed, the boost effect would cause some 
influence on the direction dependent parameters. As the earth moving velocity is 
much smaller than the light speed, the boost effect, which is about 10^'' for orbit and 
10~^ for rotation respectively, leads to little impact to the numerical values of the 
direction independent parameters {aP and c^'^) in the earth inertial frame. But the 
influence on the direction dependent parameters and , which are about 10^** 
order comparing to the corresponding direction independent parameters, needs to 
be compared with the upper bounds given in Ref. 1511 In our model 2, the direction 
dependent parameter induced by boost effect has the same order comparing 
with the upper bound in Ref. 1511 Further constraints from future experiments for 
direction dependent parameters would be a challenge for this model. Therefore 
Ref. |5T] has put a strict constraint for the sun-centered directional independent 
models built in this work. However, a model with smaller direction independent 
parameters, such as our model 1 and model 3, is consistent with the experimental 
data. 

In addition, LSND data are not used during the calculation of the three different 
kinds of models. Substituting the LV parameters numerical results of these models 
into the equation P{v^ — > Ve) respectively, we find that the oscillation probability 
are much smaller than the oscillation probability of (0.26±0.067±0.045)% reported 
by LSND. So in our work the explanation of the LSND data based on the two flavor 
neutrino oscillation is unfavored. While in this work we assume that only a limited 
special number of LV parameters are nonzero, there are still large freedoms to fit 
experimental results with Lorentz violation models, e.g., Ref. 36 explained all class 
of experimental data, including LSND, with a massive Lorentz violation model. How 
to explain LSND data would need more experimental and theoretical work. 
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4. Conclusion 

In this paper, we calculated Lorentz violation (LV) contribution to neutrino oscil- 
lation by the effective field theory for LV (Standard Model Extension). We assume 
that neutrinos arc massless and that there are only three generations of left-hand 
neutrinos and right-handed antineutrinos in nature. Unlike the conventional mas- 
sive neutrino scenario, in our calculation neutrino flavor states are not the mixing 
states of neutrino mass eigenstates but neutrino energy eigenstates. In this work, 
the equation of neutrino oscillation probabilities given by Eq. (fTT|) is a complete an- 
alytic equation and contains hundreds of LV parameters. In section 3, we introduced 
three direction-independent specific models. Comparing with neutrino experimental 
results, the upper bounds on LV parameters can be figured out. We have checked 
our calculation results with the experimental results in Ref. [511 In model 2, the 
directional dependent parameters c^^ induced by the boost effect are smaller than 
the upper bounds given in Ref. [51] but with the same order. And slightly media- 
tions for the directional independent parameters in model 2 may lead to its failure. 
So Ref. 'Wr has put a strict constraints for the sun-centered directional independent 
models built in this work. And a model with smaller directional independent param- 
eters, such as model 1 and model 3, is much better to accord with the experimental 
data. From Figs. [TJ-[5]we see that neutrinoes do not have drastic oscillations at low 
energy and oscillations still exist at high energy. Furthermore, the oscillation ampli- 
tude varies with the neutrino energy. All these characteristics are different from the 
massive neutrino scenario. In addition, the three models in this paper have different 
characteristics, too. The main difference lays in the relationship between oscillation 
probabilities and neutrino energy. In models 1 and 3, the oscillation amplitude varies 
in different energy scale and goes to zero when the neutrino energy is high enough. 
But in model 2, the oscillation amplitude is steady at high energy. These different 
models can be tested and constrained in high energy neutrino experiments. 

We also aware the neg ative rep ort for Lorentz violation from existing atmo- 
spheric neutrino data. Refs.^^^ESlM] presented detailed analysis of the zenith angle 
distribution of atmospheric neutrino events from Super-Kamiokande underground 
experiment. The analysis of super-Kamiokande data disfavors Lorentz-invariance 
as the leading source of atmospheric neutrino oscillation. There are several poten- 
tial reasons to explain this difficulty. First, Lorentz violation is not the leading 
mechanism for atmosphere neutrinos and the contribution of massive mechanism to 
neutrino oscillation can not be neglected at this range. From another point of view, 
in our calculation the mixing part is a function of neutrino energy and as Fig. [2] 
shown, the oscillation amplitudes are suppresse d at high energy, which is different 
from the Lorentz violation model used in Refs.^^^ .53i54!_ (''pj^g effect of oscillation 
suppressing at high neutrino is general consistent with the conventional massive 
model, though the suppressed threshold energy is different between the two kinds 
of mechanisms). So these novel effects in our model would produce some new effect 
in analyzing atmosphere neutrino from the view of Lorentz violation. 
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In the above discussion, we assume that neutrinos are massless and that Lorentz 
violation is the only origin ol neutrino oscillation. So we can clearly see what kind 
of new effect will occur if Lorentz violation exists in neutrino sector. It is possible 
that neutrinos have small mass and both Lorentz violation and the conventional 
oscillation mechanism contribute to neutrino oscillation. Then LV parameter will 
be further constrained and our numerical calculation for LV parameters could only 
be considered as a upper bound for LV. It is difficult to distinguish between these 
two kinds of mechanism at low energy. But high energy experiment is a good avenue 
to test Lorentz violation models. In the conventional massive model, the mixing part 
is independent of energy and neutrino oscillations disappear at high energy because 
the oscillation P is proportional to L/E. But in Lorentz violation models, mixing 
part is the function of energy and the oscillation amplitude varies with neutrino 
energy. Thus neutrino experiments with energy dependence may distinguish between 
the conventional massive neutrino scenario and the Lorentz violation scenario for 
neutrino oscillations. 
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